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Abstract 

We study the hydrodynamics of the A-B interface with finite curvature. The 
interface tension is shown to enhance both the transition velocity and the amplitudes 
of second sound. In addition, the magnetic signals emitted by the growing bubble 
are calculated, and the interaction between many growing bubbles is considered. 
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1 Introduction 



The interface between the A- and B-phase of 3 He is the only known boundary between two 
superfTuid systems. The interface motion that accompanies the first order phase transition 
from the undercooled A- into the B-phase is subject to extended experimental [1 - 4] and 
theoretical [5 - 9] investigation. 

Within the hydrodynamic regime, the structure of the interface dynamics can be rig- 
orously determined by connecting the hydrodynamic modes on both sides of the interface 
with the appropriate connecting conditions (CoCos). This has already been done [0, f| |§ 
for the planar interface. In the present work, we generalize the hydrodynamic description 
to interfaces of finite curvature, taking as an example spherical B-phase bubbles growing 
in a bath of undercooled A-phase. 

To gain a qualitative understanding, consider first planar growth. For small growth 
velocities, the released latent heat serves as the source for two second sound step functions, 
emitted in opposite directions 0. The steps move with c 2 a and — c 2 b if the A-phase is 
located on the right and c 2 denotes the second sound velocity. If the interface velocity u 
becomes larger than c 2 , both steps are left behind in the B-phase. The same is true for 
spin waves if an external magnetic field produces excess magnetization at the interface || . 
The reference velocity is now of course the spin wave velocity c s , rather than c 2 . 

Qualitatively, bubble growth is not much different, except that the surface tension 
delivers an extra push towards the energetically more favorable B-phase: Neglecting the 
anisotropy, radial step functions of second sound are emitted out- and inward for u <C c 2 . 
The inward moving step functions quickly equalize the temperature distribution within 
the B-phase bubble. For u ^> c 2 , both radial waves of second sound remain inside the 
bubble while the temperature remains homogeneous in the A-phase. The spin wave and 
the magnetization behave accordingly. (Though spin waves are emitted only for higher 
interface velocities it^ lm/s ^> c 2 |§) The case it 3> c s may seem academic but one really 
cannot discount the possibility that it was indeed already observed [[|, |8|]. Then, even both 
spin waves remain inside the bubble where neither the temperature nor the magnetization 
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is homogeneous. Causality sees to it that the A-phase is completely unperturbed. 

The complication of bubble growth is in its time dependence. While planar growth keeps 
both u and the amplitudes of the emitted waves constant, neither are in bubble growth. 
The growth velocity is not constant because the pushing surface tension, proportional to 
the inverse of the bubble radius R(t), is not. The amplitudes are not constant, because 
first of all, they vary as 1/r when moving away from the interface; second, the amplitudes 
in addition are functions of R(t). (In fact, there is a retardation effect, the amplitudes are 
functions of R(t — tf), where tf is the time of flight tf = (r — R)/c 2 or (r — R)/c s .) 

The anisotropy of c 2 and c s will distort the radial wave into an elliptical shape, with 
the principal axes being parallel and perpendicular to £, the orbital preferred direction. 
(We presume £ to be uniform.) In addition, the growth and Kapitza coefficients are also in- 
dependent, and this directly distorts the bubble. (Given enough time, of course, the surface 
tension will drive the bubble into the shape that minimizes the surface energy; generally, 
it is somewhere inbetween, reflecting some stationary balance between the surface tension 
and the uneven growth.) Assuming the same eccentricity for both ellipses, we can rescale 
to obtain isotropic growth and radial waves to perform the calculations. Afterwards, the 
distortion into elliptical waves is achieved by taking the correct value of the anisotropic 
velocity for the two principle axes. 

If there are many bubbles present, their growth can be considered completely indepen- 
dent for R > c 2 , because the magnetic signals barely feed back to the growth M, while 
the second sound waves, which do feed back, are confined to the interior of the bubbles. 
For R < c 2 however, the many emitted second sound waves heat up the surrounding A- 
phase. Since the latent heat is not completely evacuated, the growth velocity slows down 
accordingly. 

A B-phase bubble is thermodynamically stable only if its radius exceeds the critical 
value R c = a/(p A/i) R |TT|, where a denotes the interface tension and A/i = \ib — 
is the difference in the chemical potential across the interface (for every quantity X, we 
define: (X) = (X4 + Xb)/2 and AX = Xb — X&). Since A/z vanishes at the coexistence 
point, R c diverges if T approaches Tab- At lower temperatures, R c < Iqp, where Iqp is 
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the mean free path of the quasi particles, the hydrodynamic description becomes valid 
only after the bubble has nucleated and reached an appropriate size by ballistic processes. 
Fig. 1 shows the respective length scales as a function of temperature. As pointed out by 
Leggett and Yip || the nucleation process itself is likely to be of exotic nature. 

The CoCos 0, || || , we shall employ fall into three categories, those for the conserved 
quantities, those for the symmetry variables and the parameterization of the surface entropy 
production. For the locally conserved quantities, the CoCos state the continuity of the 
respective currents, expressing the fact that conserved quantities cannot accumulate within, 
or be depleted from, a region of microscopic width. This is not true, however, if the 
momentum current traverses a curved interface, where the radial bulk momentum current 
differs by the interfacial momentum current [TT| a/R. (Note that a < 0: In equilibrium 
the pressure is larger on the A-phase side.) With prime denoting the rest frame of the 
interface, the CoCos are: 

Ag' = 0, (1) 

A(p + no = |, (2) 

AQ' = 0, (3) 
Aj spin = 0. (4) 

The notations: mass current g', pressure p, nonlinear part of the momentum current II', 
energy current Q' and spin current j spin . 

The behavior of the symmetry variables in the interface region is related to the elasticity 
of the order parameter within the interface. Since the elastic coefficients are of comparable 
order of magnitude in the interface as in the bulk liquid, we have coherence of the symmetry 
variables across the interface. More specifically, we have for the phase angle <j) an d the 
rotation angle G in spin space 

A0 = 0, (5) 
A0 = 0. (6) 

It is by means of Eq. (El) that the interface tension acts in accelerating the growth process. 



Since ([5[) implies A/i = 0, the interface tension couples to the dynamics through the 
pressure dependency of the chemical potential and CoCo (@) is implicitly incorporated. 
When the growth starts at R = R c we find the initial interface velocity to be twice as fast 
as in the planar case. This effect should be accessible to experimental probe. 

The irreversibility of the phase transition is accounted for by the growth coefficient 
K for u ^> C2, and by the Kapitza resistance k for ?i <C C2, respectively. Both K and k 
are effective, integrated quantities and functions of various surface Onsager coefficients, 
thermodynamic susceptibilities and bulk transport coefficients ||. The rate of interface 
entropy production R$ = —(T) A f is then respectively given as 

a' 2 (f) 2 
Rs = 9 1 r or # 5 = — , (7) 

K K 

where /' is the entropy current. So except Eq. (|2]), all CoCos retain their form of planar 
growth. 

The solution of these CoCos in Sec. |3] shows the following result: The growth velocity 
decreases with 1/R towards the limiting value Roo (of the planar growth); while the ampli- 
tudes of second sound generated during the growth show some intricate time dependency. 

In Sec. 0, we discuss the magnetic signals emitted by the growing bubble. These are 
more accessible to experimental probe. 

The general case of N interacting bubbles V is discussed in Sec. 

2 Spin waves 

We consider the magnetic signals emitted during bubble growth for R <C c s . As discussed 
in Ref. the magnetization is for R^ (H /kOe) -1 m/s essentially given by its equilibrium 
values Xa,b H. The small deviations from the equilibrium values do not change much, going 
from planar to bubble growth. For higher velocities, R^lm/s, the dipole-interaction is 
negligible, the magnetic signals turn into spin waves and become delocalized. Here, the 
circumstances of bubble growth are more complicated and are what we shall consider. 
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2.1 Interface dynamics 

The feed-back of the spin-dynamics to the growth velocity R is negligible if R is not too 
close to c s a ||- Hence, we may take R, as calculated in Sec. [3] without the inclusion of the 
magnetic fields, as an input: 

R= Roo (i+^y ( 8 ) 

where R^ denotes the limiting value of the interface velocity for large R, i.e. R^ = it of 
planar growth. The explicit values for R^ are given in Eqs. (|27| ) and fl39|) for weak and 
strong undercooling, respectively. The solution of the differential equation (H) subject to 
the initial condition R(t = 0) = R then reads 

t Roc = R - Ro + R c In ( : |^| £ ) • (9) 
The inverse of @ can also be written in a closed form, see App. |A|. 



2.2 Spin-Hydrodynamics 

We shall work with the longitudinal model ||, which contains only two dynamical variables: 
the magnetization and the dipole angle (where 0yt is that between d and £, and 0# is 
the angle around n). 

The initial scenario contains a B-phase bubble of radius Ro, exposed to a homogeneous 
magnetic field H . Within the bubble, we assume a uniform n-texture, parallel to H; in 
the surrounding A-phase we have d _L H, i || d. Although the uniformity of the n-texture 
contradicts the equilibrium A-B interface ||12||, we shall neglect this complication as in Ref. 



The coordinates are chosen such that H points in the z-direction, while the x-axis co- 
incides with i. The spin-hydrodynamics in the A-phase is then given by the wave equation 

e A - c s p 2 x e A - c s 2 ± (d 2 y + d 2 z ) e A = o, (10) 

where c 8 » and c s± denote the spin wave velocities parallel and perpendicular to the preferred 
direction £. The consequences of neglecting spin diffusion will be discussed at the end of 
this section. 
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Eq. [10] shows that the outgoing wave fronts will be of ellipsoidal shape with a rotational 
symmetry around £. Assuming a growing bubble with the same eccentricity, as discussed 
in the introduction, we may distort the coordinates to render both ellipsoids (interface and 
wave front) spherical: 

9 a - c s 2 A A r Q A = 0, lu a = 9a, v a = <9 r 9A- (11) 

(t>A is the spin-superffuid velocity) 

The B-phase spindynamics need not be considered explicitly. After the initial time lag 
of Rq/c sB , the imploding wave front meets at the center. There, the excess magnetiza- 
tion (always parallel to the external field H) from opposite directions add up, while the 
spin superfluid velocities V0b (all having the same sign along the radial direction) from 
opposite directions cancel. So homogeneous magnetization and vanishing V 9# is a good 
approximation. 

2.3 Connecting Conditions 

With R <C c s a, we can safely neglect terms of order Rva and R5S; the CoCos (|J) and (|) 
therefore take the simplified form 

= Atu, (12) 
= jHRA X -Xac s 2 aV A , (13) 

with gyromagnetic ratio 7 and susceptibility \- Working with the longitudinal model, only 
one of the three relations (||) is needed here. 

2.4 Solution of the Connecting Conditions 

The procedure to determine the space and time dependency of the radial waves from their 
CoCos is given below: 

In the first step, we calculate the values of the general solution 

9a = /(r ~ Cs ^ } (14) 



at the interface (subscript I) from the CoCos ([L^, |13|). 

Substituting the time dependent bubble radius R(t) for r, we find 

f'(r-c sA t) 



ua\i = (d t e A (r,t)) \ r=R{t) = -c sA 

ff'(r-c BA t) f(r-c sA t) 



(15) 

r=R(t) 



v A \j = (d r Q A (r,t)) \ r=R(t) 



(16) 

r=R(t) 



where /' denotes the derivative of / with respect to its argument. Now, the amplitudes 
at the interface are functions of time only. For sake of clarity, we introduce an auxiliary 
function by 

F(t) = f(R(t)-c BA t) = f(R{t)), (17) 

implying 

f-g®-. (18) 

K — C sA 

The CoCos ( |T2| , |T3| ) thereby turn into differential equations for F(t) which can be solved. 

Having obtained F(t) with the appropriate initial condition F(t = 0) = 0, the next step 
is to reconstruct the (r, t)-dependency of /, where we need the inverse function of R(t): 

F(t(R))=f(R). (19) 



Given f(R), we can substitute R = r — c sA t and obtain u A and v A for all r, t via Eq. (|TT|). 

For technical reasons, however, we modify this procedure and eliminate the time t with 
the inverse of R(t), instead of eliminating r by R(t), as described above. The reader is 
referred to App. ^ for the details of the calculation and the result. 

For the qualitative features cf. Figs. 2 and 3. The spatial variation of the outgoing 
radial wave is illustrated for different times. Each curve displays two effects. First the 
1/r-decrease as the wave moves out, and second the amplitude at r, t is a function of R 
at t — (r — R)/c, a retardation effect. The line at the upper surface of the box marks 
the location of the bubble surface, relative to the starting value. Below, one can track the 
decrease of the amplitudes at the interface. As will be shown in Sec. |3|, these plots also 
visualize the second sound amplitudes, emitted by a slowly growing bubble. 
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At the head of the shock wave, i.e. where r = Rq + c sA t, the expressions for the 
amplitudes simplify to become 

v A {t) = 5— - , u A {t) = -c sA v A (t). (20) 

c sA X c sA -"-0 + C S A 1 

The qualitative effect of the spin diffusion, neglected here, is the rounding of the mag- 
netization step with out at the wave front. So Eq. (p0| ) is valid for the total change in va 
and co a. 



3 Growth Dynamics of spherical bubbles 

The interface dynamics of superfluid 3 He is fundamentally different from other systems 
jTj. If u <C c 2 , the evacuation of latent heat is very efficiently accomplished via second 
sound, the difference between hypercooling and normal supercooling is inessential, and u 
is determined by the Kapitza resistance of the interface. For ii 3> c 2 , the system is always 
hypercooled, no heat evacuation needs to take place and not much does, u is determined 
by the growth coefficient. 

In bubble growth, two similar scenarios prevail: For weak undercooling of the A-phase, 
R <C C2, an outgoing radial wave is emitted into the A-phase, while an imploding wave 
is emitted into the bubble. The latter will interfere at the center of the bubble, as in the 
case of spin waves, yielding vanishing counterflow and an averaged temperature, depending 
only on the current bubble radius. 

For strong undercooling, assuring R 3> c 2 at all time, no sound is emitted into the 
A-phase, while an exploding and an imploding wave superpose inside the bubble, yielding 
a time dependent temperature distribution. 



3.1 Bulk Hydrodynamics 

We start from the linearized equations of bulk hydrodynamics (correct for both limits 
I > c 2 , i? < c 2 but not in-between, R w c 2 [0), in a form appropriate for the 
spherical symmetry of the problem. As in planar growth, first sound is too fast to need 



explicit consideration. Its only role is to maintain uniform pressure on both sides of the 
interface, while keeping the discontinuity Ap = a/ R across it. This leaves two equations: 
Conservation of entropy and the equation for the counterflow. The counterflow, w = w(r) e r 
is curl-free and can be written as w(r) = d r <p. The two hydrodynamic equations then 
become 

^5T + clA r 4> = 0, (21) 
d r l^+ P -^5T^ = 0. (22) 

Eq. fl2"2] ) seems to allow <fi and —pa/p"5T differing by an explicit function of time. This 
is ruled out by the wave equation that can be obtained from the combination of (^1]) and 



(f2~2|), which must not have any time dependent source. 
In summary, we end up with 

n n . 

w = d r <j)(r,t), 5T = -—<p(r,t), 0-c^A r = O, (23) 

pa 

in close analogy to flTTp. 



3.2 Weak undercooling 

We will consider two different stages: The first stage is characterized by the existence of 
an incoming and an outgoing second sound wave and lasts only for a time span of order 
R0/C2, where Ro = R(t = 0). The second stage starts after the interference has established 
homogeneity inside the bubble, making the temperature solely dependent of the bubble 
radius. In the calculation of the second stage we shall neglect the existence of the first. This 
only changes the initial condition and hence the exact shape of the radial wave front but 
not its magnitude. For very weak undercooling, R c is large compared to the quasiparticles 
mean free path, so that for t = R = R c and T = T« are the valid initial conditions. 
Otherwise, we assume a nonvanishing initial temperature discontinuity A Tj = T^b — TiA 
and R = R > R c at t = 0. 
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We expand the CoCos (|3],|5],|7|) , t° first order in R, 5Ta,b, wa,b around the initial tem- 
perature Tf. 



= -pRAtJi + p s (a) (w B -w a ), 

= Ap t ,(l + ^)-(a}(5T B + AT l -5T A ), 







-pR(a) - k{5T b + ATj - 6T A ) 



(24) 
(25) 
(26) 



all quantities refer to the laboratory frame. Terms of Order (A<Ji/{a)) 2 as well as d p (p s / p), 
dr(p s / p) and Ap/p are omitted. 

From (|25|) and (^) we readily obtain the interface motion as quoted in Eq. The 
limiting value of the interface velocity is 

k Api 



Ra 



(a)*p 



(27) 



We consider the sound amplitudes at the interface for the first stage. Starting with 

,/ ,x f(r-c 2 t) gjr-cyt) 
<P{r,t) = 1 , 

we apply the method described in Sec. ^| and calculate the sound amplitudes to lowest 
order in (R — R c )/R c : 



5T 



A/ij 1 . m Roc A(Ti 
zr LZ _|_ -AT- — — 

T (a) + 2^ U c 2 (8 T a) 

4.1 S ^ti ( ° 2 4- 1^ 1 at ° 2 - ;Roo^;\ R-R c 

2 I (*) Uoo 2 ^ c 2 (a) J i? c 



(28) 



i? — R r \ 2 



B 



p'"C 2 

1 [A// 



(a) + 2 " + c 2 («9 T a> 



(29) 



+ 

+0 



2 |» V /? 

i? — R, 



±2-^ + 3 — ATj (±-¥- + 1 



R r 



oc 

2\ 1 



R„ 



Aoj 
(d T cr) 



l±3- 



Rr 



C 2 



R — R c 

Rr 



Setting ATj and R — R c to zero, these expressions can be compared to the results of 
planar growth 0: As it turns out, the surface tension doubles the initially emitted sound 
amplitudes. Note that the terms of first order in R/R c —1 are dominated by the contribution 
proportional to CijRoo- 
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For the second stage, the situation is quite similar to that of spinwaves, examined in 
the last section. Since (|25f) and ( pjf) have already been used once to determine the bubble 
radius as a function of time, we shall now only consider the CoCos (^) and (|25|). With 
Wb vanishing, (|24D simplifies to 

= pAciR + p"(cr)w A , (30) 



which is analog to (p~3|) when wa is identified with —v A - Combining ( |30"D with fl2~5|) yields 



^ = TXpS ( ST b - + AT,) . (31) 



By identification of with u>a, (plQ is analogous to dl2|) . (These identifications are justi- 
fied by inspection of the bulk hydrodynamics: (p3|) goes over into ( |TTD by the substitution 
of — Qa for 0; but this does not alter the solution of the wave equation.) 
We can therefore take over the results of Sec. by setting 

w A = -v A , (32) 

ST A = -y-r — UJ A , (33) 
PS?) c sA 

when we in addition observe the respective redefinition of A and ca, 

Acnp 8 - C 2A 
A = — i— — -Roo R c , c A = -r— , (34) 

W P Roo 

introduced in App. |B|. As a consequence of this analogy, Figs. 2 and 3 also serve as 
illustrations of counterflow and temperature amplitudes. 



3.3 Strong undercooling 

At the strong degree of undercooling we are considering here, the hydrodynamic description 
only holds after the bubble has grown to an appropriate multiple of the critical size. We 
therefore impose the initial condition R(t = 0) = XqR c with Xq ^> 1. 

Inside the bubble, an incoming and an outgoing second sound wave superpose. After 
the incoming wave front has been "reflected" at the center of the bubble, there are two 
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spherical discontinuities of temperature and counterflow, expanding outward and laging 
behind the third discontinuity, at the bubble surface. 

Neglecting terms of the order wb/R, an expansion of (0,^,0), again in R c /R — 1, yields 
for the discontinuities at the bubble surface 

= TAvi + A/i, (l + -£) + Ti d T a B (ST B + AT,) , (35) 



R 



(R \ o n 
1 + ~ °b (5T B + ATJ + P — Rw B , (36) 

,,H TiAoi + (Wtcjb + ^P) (6T B + ATi) . (37) 



K V 2 

From (|35|) we obtain the temperature 



which, together with (p7|), again leads to (||). Here the limiting value of the interface 
velocity is 

By calculation of the total counterflow amplitude 

p jAfii ( a B \ a B Aai 1 1 

our initial neglect of 0{w B /R) proves to be appropriate. 



4 Interacting Bubbles 

As mentioned in the introduction, we expect the growth velocity of the bubble to slow down, 
when its surrounding heats up from second sound radiation sent out by other bubbles. Only 
the case of weak undercooling (R < cq) is worth studying, since for strong undercooling 
(R > C2) all the latent heat remains inside the bubble. 

Let us for simplicity consider a volume V of A-phase, containing iV B-phase bubbles, 
nucleated simultaneously and distributed uniformly over V. Each individual bubble grows 
as described in the preceding section until it is reached by the second sound pulses of its 
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neighbors. From now on, the accumulative heating of the A-phase increasingly hinders the 
evacuation of latent heat and slows down the growth velocity. 

We will restrict our calculations to the regime where the temperatures of both the 
A- and B-phase can be regarded as being uniform. The Volume V is thus divided up 
into two subsystems: iV bubbles of total volume V Q = N4tc/3R 3 and the remaining A- 
phase of volume V — V D . Since the gradients vanish, the local hydrodynamic description of 
the previous sections reduce to thermodynamic considerations: For small deviations from 
equilibrium, energy conservation accounts for equality of the entropy currents, 

S B = ~S A , (41) 

where Sb = V Q pas and Sa = (V — V ) per a denote the total entropies of the respective 
phase. Because of the identity 

d t S B = jdV & B + Jds Rob = - jds f B , (42) 

Vo dVo dVo 



we identify Eq. (f|I]) as the continuity condition of the local entropy current, integrated 
over the area of the interface. In the same manner, we obtain the Onsager relation, starting 
from f' B = kAT as 

-^f = n(ST B -5T A ). (43) 
At 

Here, Aj denotes the area of the interface so that the total rate of interface entropy pro- 
duction is an extensive function with respect to the interface ctrGcL, clS it must be. Being of 
intensive nature, the continuity condition for the phase angle 

= Afa (l + ^ ) " °b ST B + a A ST A (44) 

does not change its appearance. 

Because of the homogeneity of the subsystems, all thermodynamic quantities are explicit 
functions of R. We therefore need to specify only a single initial condition: 

6T B (Ri) - 6T A {Ri) = AT,, (45) 
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where Ri denotes the initial bubble radius and A T« is the temperature difference that was 
built up across the interface during the switch-on process. From Sec. [3|, we expect A T, L 
to be of the order A/Xj/(er). Since (|4]) was already valid before homogeneity had been 
established, the initial temperature amplitudes are found to be 

*r A m = A »-'"" r : sr„W = *"-° A * T ' . (46) 

Here and below, we neglect the effect of interface tension, i.e. we assume R c <C Ri. 

From ([II]) and the derivative of (|44| ) with respect to R we obtain 3r5Ta and dn5T B 
which may be integrated to give the temperature amplitudes as functions of the bubble 
radius: 

tfZ^) = STab {Ri) - — ^— In ( 4ff/3M 3 fl + ]/& J . (47) 

with a = cr^ c^erg — o"b ^t^a, b = a B Qt&a- 
The growth velocity follows from 



p tr B + Ot(?b R d R 5T B /3 

Obviously, the system approaches an equilibrium state with equal temperatures on both 
sides of the interface and vanishing R. Starting with 

the interface velocity decreases and eventually comes to rest at the final bubble radius Rf. 
The B-phase bubbles then occupy the volume 

and the final temperature amplitudes amount to 

ST A (R f ) = ST B (Rf) = (50) 
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A Calculation of x(r) 

For convenience, we introduce reduced dimensionless variables by x = R/R c and r = t/t c , 
where t c = R c /Roo is the characteristic time for interface movement. The solution (Qj of 
the equation of motion of the bubble radius (§) then takes the form 

r W = ,- Io + l n fi±^V (51) 



1 + x 

which shall now be inverted to yield the bubble radius as a function of time. 



The inverse of z i— > z e 2 , is known as Lambert's function W BM- For real valued x, W 



has two branches: 

W (x) > -1 ; -- < x 
e 

W-i(x) < -1 ; — < x < 
e 

Exponentiation of (|51| ) yields 

-(1 +x)e~ {l+x) = -(l+x )e- (1+xo) e- T . 

By definition of W this is equivalent to 

x(t) = -1 - (-(1 + x )e- {1+Xo) e- T ) 

W-i is the right branch, because the argument is always negative, and x must not be 
limited. 

Note that we reserve x and r to denote the reduced bubble radius and the corresponding 
reduced time, respectively. Arbitrary coordinates (r, t) are referred to by (y,$), see App. 
B. 

B Spin wave amplitudes 

With y = r/Rc, ■& = t/t c , <p(y,&) = 9 A (r,t), v A = v A R c , u A = uj a R c /c sA , c A = c sA / R^, 
Eqs. fln|) read 

v A = d y if, Uj a = d# if, dlf - c 2 A A y if = 0, (52) 
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while the CoCos ( |T2| ) and ([13]) become 

= Atu, = A d T x + v A , with A 



-Roo R r > 0. 



(53 



The general solution of (f)%) is </>(?/, ■#) = F(y — c A $)/y- Inserting v A \j = F' \ T /x — F \ x / 



x 



d x q{x) 



and uja \j = —F' \ j /x, with 

F \i = F(x - c A r{x)) = q(x), F' \, = " XHy "' , 

1 — CaO t X 

into the CoCos (0), we obtain the following differential equation for q(x): 
d x q{x) = X(1 ~ C ^\ (X) - A (x (1 - c A ) + 1) • 

X{± + X) 

The special solution of (ESI) subject to the initial condition q(x ) = is 



(54) 



(55) 



q(x) 
C(x) 



-A- 



x 



1 + x 



(■A 



(C(z) - C{x )) with 



(56) 



1-ca 



1 + x 



,1+CA 



+ —(1 + x) CA 2 F 1 (-c A , 1; -c A + 1; 



1 — C-A C\ '' 1 + X 

Here, 2 F\ denotes the hypergeometric function as defined in |I3J, Eq. (7-76). Choosin 
Xi(y, $) so that, 

F ( Xl (y, $)-c a t ( Xl (y, 0))) = F(y - c A 0), (57 
we obtain the spin wave amplitudes outside the bubble as 



VA{y,$) 



d x q(x) 



y 1 - c A d x r(x) 
1 d x q(x) 



x=x\(y 



y 



2 ( l\ X ) \x=x 1 (y,tf) 



y 1 - c A d x r(x) 



(58 
(59 



x=xi(y,tf) 



A closed form expression for xi(y,$) is derived in App. |C[ 



C Calculation of X\(y, $) 



A necessary condition for (|57|) to be fulfilled is given by 

xi - c A r(x 1 ) = y-c A {}. 
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Inserting the definition (|D of t(x) gives after exponentiation 

l ~ c A f . . v i^A(l+ 2x ) l-C A/ IzlA v^A± 

(1 + xi) e c a y Ll — (1 + x ) e c a e e c a , 

so that 

Xl (y, tf) = -1 + f — ^(1 + s ) e^T e"^ e^T^) . 

1 - ca V ca / 

To check for the right branch, we make use of the identity 

x = xx(x,t(x)) 

1 c A (\-c A izfA (1+x) \ 

= + - W-i\ (l + x)e c a K . 

1 - ca V c A J 

Since > 1, the argument of W-\ is negative. W_i is the right branch, because x has to 

be a monotonic function of itself. 
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Figure Captions 



Figure 1: Mean free path Iqp of quasi particles in the B-phase (taken from Fig. 10.1 
of Ref. |10[) and critical bubble radius R c (as given by Kaul and Kleinert JTT[),both as 
functions of temperature. 



Figure 2: Counterflow and spin-superfluid velocity, scaled as 



Figure 3: Temperature and excess magnetization, scaled as 
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